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Abstract This work proposes a very simple random matrix model, the Flip 
Matrix Model, liable to approximate the behavior of a two dimensional elec- 
tron in a weak random potential. Its construction is based on a phase space 
analysis, a suitable discretization and a simplification of the true model. 
The density of states of this model is investigated using the supersymmetric 
method and shown to be given, in the limit of large size of the matrix by the 
usual Wigner's semi-circle law. 

I Introduction 

This paper is devoted to the rigorous study of a random matrix model that is 
liable to give a good approximation for the density of states (DOS) of the two 
dimensional Anderson model in the weak disorder regime and in the vicinity 
of the Fermi level. The main result is that the DOS of this random matrix 
model converges to the usual Wigner's semi-circle law in the small coupling 
limit. The Anderson model of an electron in a random potential corresponds 
to the Hamiltonian 

H = (-A + XV(x)) 

acting on the Hilbert space L 2 (M 2 ), where A is the usual Laplacian and V is 
a real Gaussian process on R 2 with short range correlations. 
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This model was initially proposed to describe the electron motion in doped 
semi-conductors at low temperature or in normal disordered metals. It is now 
the central model for the theory of eletronic transport and wave propagation 
in disordered systems [p]]. It was conjectured by Anderson as soon as 1958 
|2| that such a model exhibits a localized phase in which the electrons are 
trapped by the defects. In 1979 it was argued that this model has a phase 
transition in dimensions three or more between the localized phase and an 
extended one ||. 

The localized phase is now well under control. In one dimension, localization 
was rigorously established for any disorder at the end of the seventies j|, [5| . 
Later localization was established in any dimension at strong disorder or for 
energies out of the conduction bands ||. A simplified and more efficient 
method to get this result was given in J7|]. 

In contrast the weak disorder regime is still poorly understood. In two di- 
mensions it has been argued || [8| and numerically established || that lo- 
calization persists at arbitrarily small disorder, with a localization length 
diverging like 0(e C//A ). In dimension three, numerical simulations confirm 
the existence of the Anderson transition |J, leading to an extended phase. In 
addition, analytical results |10| and other numerical calculations show that 
the level spacing distribution follows the Wigner-Dyson distribution for ran- 
dom matrix theory (RMT) (llj. This gave the motivation for a description 
of mesoscopic system in terms of RMT [[HJ. This method has been very suc- 
cessful when compared to experiments and was the source of developments of 
super symmetric methods []T3|, [HJ in solid state physics. But on the rigorous 
level it is still a mathematical challenge up to now to prove even regularity of 
the DOS (e.g. real analyticity in energy in a certain interval) at arbitrarily 
weak coupling in dimensions two and three. 

Using a phase-space analysis inspired by the renormalization group method 



around Fermi surfaces in condensed matter |15|, 16, I7| , G. Poirot and coau 



thors [|18| , [19f| have established that the effective Hamiltonian near the Fermi 
level is given indeed by a random matrix model. In two dimensions this ran- 
dom matrix model is similar to the GUE, but contains an extra discrete sym- 
metry. Unfortunately in three dimensions and beyond, this extra symmetry 
becomes continuous, and produces more complicated correlations between 
matrix elements [0. Recently a band matrix model in three dimensions 
has been proposed and its DOS studied through the supersymmetric method 
1 201 . In this note, motivated by these works, the two dimensional case for 
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the Anderson model is reconsidered. A simplified model is constructed on 
the basis of the phase space analysis. It is a matrix model, although not 
one of the regular ensembles. We analyze this model rigorously through the 
super symmetric method, and proves that its DOS converges to the Wigner 
semi-circle law as the matrix gets large. 



II The Flip Matrix Model 

Average spectral properties can be studied through the averaged Green's 
functions of the model. For some suitable choice of units, these averaged 
Green's functions are defined as 

G±(y-x) = lim [ du(V)(x\(-A-E + ie + XVyHy) 

6^0± J 

where dfi(V) is the disorder distribution, A is a coupling constant control- 
ling the disorder strength, and E is complex with small imaginary part e. 
The DOS is defined by the imaginary part of the averaged retarded Green's 
function : 

v(E) = -ImG+(0) 

7T 

For a given non-zero E in the conduction band (e.g. for \E\ neither too 
large nor too small), the free Green's function (p 2 — E)~ l (we simply write 
—A = p 2 ) is singular on a surface, which in two dimensions is a curve. Two 
regimes have been clearly identified: 

1. for \p 2 — E\ > 0(1) A 2 , the random potential XV is indeed a weak pertur- 
bation of the free Green function. It is statistically unlikely to develop 
an eigenvalue of the same size than p 2 — E, making the denominator of 
the Green's function singular. So this regime can be controlled by the 
usual techniques of multiscale cluster and Mayer expansion, with some 
large/small field analysis. Typical V's belong to the small field regions; 
exceptional Vs in large field regions are controlled through some type 
of Tchebycheff inequalities, and give rise to small probabilistic factors. 
These factors can then pay for the necessary complex contour defor- 
mations and rough bounds that desingularize the denominator of the 
Green's function in these regions |T9 | . 
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2. for \p 2 — E\ < 0(1) A 2 , the interesting region, a non trivial imaginary 
part generated by the potential should stop the renormalization group 
analysis. This is easily seen at second order perturbation theory in 
A: an imaginary term proportional to iX 2 appears. Hence G± should 
decay exponentially with a length scale proportional to A~ 2 , and the 
DOS should be regular. The problem is to justify non pertubatively 
this well-known theoretical argument. 



For the moment the main result in this direction is 21 : 



Theorem 1 There exists some finite rj > such that for e = X 2+v the 
averaged Green's function < G +t (x,y) >y decays at large \x — y\ with a rate 
r(A) which is independent off], and inverse of X 2 . 

This result was obtained through some difficult non-perturbative Ward iden- 
tity which uses the particular momentum conservation laws in two dimen- 
sions. It is non-trivial, since the decay rate in A -2 , as expected from pertur- 
bation theory, does not depend on the 77 cutoff. This result seems difficult to 
extend to fully perform the limit e —>■ (this would in fact be implied by a 
similar theorem but with r] = 1). 

To complete the proof of exponential decay of the averaged Green's function 
and of the regularity of the averaged DOS (without X 2+v cutoff), the super- 
symmetric method [JT3], [TJ]] may be more convenient. Indeed the \p 2 — E\ < 



0(1) A 2 regime of the 2D Anderson model is no longer perturbative in the naive 
sense, since the potential combines non-trivially with the deterministic part. 
Supersymmetry seems the right tool to control this phenomenon since in this 
formalism it corresponds simply to contour deformation to a non-trivial sad- 
dle point that generates Wigner's law and stops the RG flow. The functional 
integral away from the saddle point can be controlled non-perturbatively by 
a small/large field analysis, which is compatible with standard weak-coupling 
cluster expansions p0| . 

The Flip Matrix Model studied in this paper incorporates four simplifications 
with respect to the true Anderson 2 dimensional model: 

• The regime \p 2 — E\ > 0(1)A 2 has been completely removed from the 
model. This is justified since this regime being fully perturbative can 



be added later in the style of [19|. So we restrict our Hilbert space 



to functions supported in momentum space by a tiny shell around the 
circle p 2 = E (or 2 — cospi — cosp2 = E for a square lattice version of 
the model). 
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The space has been reduced to a single cube of side size A~ 2 . This is 
in agreement with the idea of a contour translation generating a i\ 2 
part in the denominator of the Green's function. Then the full model 
should correspond to weakly coupled such cubes. Once the single cube- 
model has been understood, the full model should follow from a cluster 
expansion. 

In this cube the momentum shell has been divided into 0(A -2 ) cells, 
called sectors, and the corresponding infinite dimensional Hilbert space 
has been replaced by a finite dimensional one, so that the random 
operator corresponding to V is replaced by a random matrix. This 
discretization step is justified since the phase space analysis is now fine 
enough, so that each random operator being approximately constant, 
it is legitimate to replace it by a single random coefficient. 

The probability law for the random coefficients of the matrix has been 
computed from the 2D momentum conservation rule, with a last sim- 
plification. As well known the rhombus rule of 2D momentum con- 
servation for four vectors of equal length slightly weakens for almost 
degenerate flat rhombuses and this degeneracy requires an anisotropic 



slicing of sectors and of the cube into parallelepipeds [[u], [19], We 
forget this complication to stick with simpler isotropic sectors, so we 
use a simplified momentum conservation rule. Normally introducing 
anisotropic analysis is only a (painful) complication of the model that 
gets rid of the difficulties associated to almost degenerate rhombuses. 

Then the subsequent model would be almost the GUE if not for the rhombus 
flip. This flip creates a new symmetry in the random matrix, so that the 
supersymmetric analysis of this model requires essentially twice as many 
variables as for the GUE. The full proof of regularity of the DOS for the 2D 
Anderson model is therefore now reduced to the hard task of fitting together 
all pieces of the puzzle by removing our four approximations to go from the 
Flip Model to the true 2D Anderson Model. 



Ill The Flip Matrix Model: notations 

Let N be a large number parametrizing the number of sectors. To model 
the 2D Anderson model in a large spatial cell of side size 0(X~ 2 ) ,we have 
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divided as indicated above the momentum shell of width 0(A 2 ) around p 2 = 1 
into 2N = 0(X~ 2 ) isotropic cells. Since the initial V(x) was real, its Fourier 
transform V(p) (the hat on Fourier transforms is removed for simplicity) is 
complex and obeys 

V(-p) = V(p) . (III.l) 

Since V(x) is random independent identically distributed in direct space, 
V(p) is Gaussian white noise, so that 

(V(p)V(q)) = 5(p-q) . 

In order to implement the momentum conservation constraint at the vertex 
the "rhombus rule" will be simplified by forgetting the degeneracy of the 
rhombus near collapse. Discretization of V(p) follows from it being consid- 
ered as a random matrix between the cell indices. Collision with the initial 
potential V(p) transforms a state with momentum r into a state with mo- 
mentum r+p. Therefore the discretized matrix corresponds to V with matrix 
elements Vf^y) where 7 is some discretized version of p. Its coefficients are 

where a and (3 run over the set of sectors and V a $ transforms a state with 
momentum sector a into a state with momentum sector (3 if and only if 7 ~ 
(3 — a. Solving this equation for a known 7 gives in a generic case < I7I < 2 
two possible pairs forming a rhombus, namely (a, (3) and (—(3, —a), so that 
for (3 7^ ±a there is a single random variable for the two pairs: 

V a>l 3 = V-p,-* (III.2) 

In the degenerate case 7 = 0, it gives 2N equal pairs V = V a>a , and in the 
degenerate cases I7I = 2 it gives a single pair V a - a . If, as explained in the 
introduction, the almost degenerate rhombuses are neglected, which should 
be harmless since they are not generic, as well as the possible effect of nearest 
neighbor slight overlaps between the momentum definition of cells which may 
lead to slight complications, V is well approximated by a random Gaussian 
2N x 2N complex matrix. This matrix is Hermitian because of 

V a , p = Vp, a (III.3) 
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It is therefore convenient to label sectors as {1, • • • , N} U {— N, 
that {!,••• , N} labels projective sectors, as shown in Figure 1. 



■1}, 



so 




Figure 1. Numbering sectors 

This leads to a set of 2N(N — 1)/2 = N(N — 1) complex variables V a $ for a G 
{l,-- - ,iV},and/3e {1, • • • , JV}U{-JV, • • • , -l},a< \f3\. They fill a quarter 
of the matrix, made of the upper half of the upper-left quarter and of the 
lower half of the upper-right quarter. All other coefficients are deduced from 
these. Thanks to ( |111.2| )-( [I11.3 ), and to the particular numbering of Figure 1, 
there is symmetry around the antidiagonal, and conjugate symmetry around 
the ordinary diagonal (so that V is hermitian). 

The upper anti-diagonal is made of N additional independent complex vari- 
ables V a _ a for a G {1, • - • , N}, corresponding to momentum transfers of 
length 2 of V, while the lower antidiagonal is hermitian conjugate to the 
upper one. Finally the variables V aa are all equal to a single real variable Vq 
and that fills the ordinary diagonal of the matrix V. 

The total number of real independent variables among ther matrix elements, 
is therefore 2N(N - 1) + 2N + 1 = 1 + 2N 2 . They all have covariance 
equal to 1. The corresponding normalized Gaussian measure will be denoted 
dfi(V). This model is similar to the standard GUE (that would have 4N 2 
independent real variables for a 2iV x 2N matrix instead), except for the 
change on the diagonals and the additional symmetry with respect to the 
anti-diagonal, corresponding to rhombus flips. 

The structure of V can be summarized as follows (remember that the 
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rows and columns indices are labelled as {1, • • • ,N, —N, • • • , —1} hence not 
as {!,••• , N,-l,--- ,-N}): 



V 



V 



V 



v„ 



Vo 



Vo 



v» 



Our main result in this note is: 



Vo 



V 



Theorem 2 In the large N -limit, the DOS of this Flip Matrix Model con- 
verges to Wigner's semi-circle distribution. The corrections to the limit are 
uniformly bounded as 0(1/ N) as N —>■ oo. 

The rest of the paper is devoted to the proof of this result. Heuristically, 
the diagonals should not matter much for the average spectral properties. 
The main difficulty is to check that the additional symmetry of the rhombus 
flip does not alter Wigner's semi-circle law. This is made precise using the 



super symmetric formalism [13], [14 . 



IV Mean field analysis 

Let the superfields be 



^ act ( S{j a , X(TQ ) 



for a = ±1 and greek variables such as a, (3 now running only within [1, N\. 



The conventions are taken from [[Uj]. Hence bosonic fields are complex with 
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S = S + and pairs of anticommuting variables Xi X + satisfy: 

v-v- = — V-V- V^V • = — V-V^ V^V^ = — V^V^ 
A»Aj AjAi ) Aj AJ AjAj ) At Aj Aj At 

{Xi ) Xi ) (XiXj) Xj Xj i 

J " dxi = J dxf = 0, J XidXi = J xt d X. 



"2- 



so that 



/ \{dxtdx l e- x+Mx = det(M/27r) 
whereas ordinary complex bosonic variables give 

/ Y[dS+dS ie - s+MS = det(M/2 7 r)- 1 
J i 

(and M has to have a positive real part). By the usual rules of anticommuting 
integrals each Gaussian term is written as an integral over superfields. Hence 
the density of states can be written as 

v{E) = lim [s+S a e iLo T\cM+ a cM aa dv(V) (IV.4) 

ImE~*G+ IT J 

where the supersymmetric action decomposes as 

L = {A + B + C ) 



with 



a</3,<7 

B = ^\V a ,- a (*+*- a + h.c.) 

a 

C = E + XVo^Ka^aa 
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correspond respectively to the main part, the anti-diagonal and the diago- 
nal of the matrix. The Gaussian integration over the V variables will be 
performed except for Vo leading to the quartic terms 

u(E)= lim -Im [ S+S a e- L T[ d^ a d^ aa d^ (V ) 

Q;,cr 

L = \ 2 {A + B + C) 
A=J2 + ^*-«)(*J»*« + ( IV - 5 ) 

a</3,a 

a 

The sum in ( |IV.5|) can be written as 

A = J^Kp.a^a With $ Q/3;(7 = (*+* a + ^ a *_^) . 
a 

Exchanging a and (3 in A = ^aB-a^a/3;a, both for a — ±1 does not change 
the sum. This is because the bosonic part in $ commute and the rule for the 
fermionic part leads to <£^g. + = ^/3a;+' whereas for a = —1, ^p a> - — ®a/3,-- 
Hence 

A = \ E + ^*-«)(*5»*« + *la*-/j) + ^++ + D +- + ^- » 

where 

D_ = -(1/2) + *i*-a)(*l«*« + = ~2S 

Q 

D ++ = -(1/2) ^(M/+ a M/ CTCt ) 2 D + _ = - Y,^*^-*) 

a, a a 
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D+- = - ^2(S+S a + XaXa)(S-a S -a + X-aX-a) ■ 
a 

The decomposition of -D+_ into its boson and fermion parts and adding B 
above, eventually leads to L = X 2 (A + B + C) with 

A = (1/2) J2 + * + V±a*-*0) 

a 

cr CT cr 

C = C = (A- 1 ^ + V ) J2 (IV.6) 

a,cr 

The term i3 being a sum of diagonal quartic terms, will be neglected for 
a while. For indeed it cannot be written as the square of a sum over a. 
However, the bosonic part of this term has the wrong sign, which may create 
some difficulties when performing the integration in ( IV. 4 ). This problem will 



be addressed in section VI. The other terms, however, can be reorganized so 
as to get the square of a sum over the a's by pushing the terms with index 
a on the left and the ones with index j3 on the right. This is possible thanks 
to the commutation rules for bosons and for fermions. The calculation is 
tedious but straightforward and gives: 

A = -(l/^E^JKl/^E^Kk.) (IV.7) 

Q;,cr a,cr 



an 



The squares can be unfolded by mean of an integration over auxiliary gaussian 
fields. This amount to introduce two real fields ao and bo, one complex field 
a, a and two pairs of fermionic fields £,£ and rj, fj with Gaussian measure 

dfx(a , b , a, a, £ £, fj, rj) = e -^/2-^/2-J8/2-[a[V2-€'€-^^ x . . . 

dV da db d 2 a * , . . . 

■ ■ ■ ( 2 7t) 5 / 2 V V ' ^ ^ 
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Therefore 



with 



XaaXcra 



e \ 2 (A+C) = J dfJl e iAE a *±R* a 

= (a + V + \- 1 E)J2sLS™ 

a a,cr 

+ (ibo + V + X^E) 

a,cr 

+ S a S- a + h.c^j 

+ I r E S ™X*a + h.C. J + I rf J] SaaX-aa + h.C. 
\ era / \ era 

where $ is the superfield = (S£, S- a , Xa , i>X-a) and R is a 4 x 4 super- 
matrix. Setting A = (a + V + \~ l E) and iB = (tb + V + A" 1 -/?) this 
matrix is given by 



(IV.9) 



R 



( A a e -vn \ 

a A rj* -ill 

i r] iB 

\ -irf -<* -z.B / 



A p* 
p B 



(IV. 10) 



In order that the integration over the primitive fields be given in term of $ 
it will be convenient to introduce the new fermionic fields 



X—a — l X-a i X—a ~ l X-a , d\-a — ^X-a i ^X-a ~ ^X-a 



This leads to the formula 



J] / d 2 S a d 2 S^d X + a d X adxt a dx-ae lX ^ m - = [Sdet(t\R)]- N . 

a J 

The super determinant can be computed from the Schur complement formula 
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Sdet { A p b) = dib^-"*^ 1 ")- 

Defining a = a + V + A _1 -E, 6 — + Vo + A _1 i? and using the Pauli 
matrices <7 2 , cr 2 , cr 3 the superdeterminant of R is given by gives 



SdetR = — — det ( a l + a\0\ + 0202 (IV. 11) 

' c -it] \f 1 \ / e 17 \ t-i 
77* ; v 1 / v -<* J 

- {a + 2r^o 1 }{« + 2^ 1 

Including back the scaling factors, taking into account the relation 2 A 2 = N^ 1 
and ignoring the B term, the density of states is given by 



/ 



da dV ■ ■ ■ ( e 



-(al+vt+Wf+bl) (E + Vo + ib ) 2 



(E + V + a ) 2 - |a 

x(l + O(l/A0) (F ^l° + "° 1 ,2 • (IV12) 
(E + V + a ) 2 - \a\ 2 



For indeed, taking into account the part of the 2-point function involved in 
the integral ( [IV. 4| ) , the fermionic part gives a 0(1/ N) correction in the form 



N-ls ■ N + 1 



((l + rr}*rjCO ) = 1 ~ r 



AN 2 

thanks to the integral 



7T 



d£*dZe- 2NV *CZ 



2N / * * * ' AN 
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V The saddle points at leading order 



In this section the 0(1/N) and B terms are still neglected. As iV — > oo, the 
integral ( EVI21) can be treated by the saddle point method, since no fermion 
variable are involved anymore. The saddle point equations can be written 
dS = where S is the effective action given by 



S = -(al + V 2 + \a\ 2 + b 2 ) + 2 ln(E + V + ib Q ) - \n{{E + V + a ) 2 - \a\ 2 )) . 

More precisely, differentiating with respect to the independant variables ap- 
pearing in the integral ( [IV.12|) gives 



-(1/2) 



dS 
da 



a + 



a + V + E 



(E + V + a ) 2 - \a\ 2 



"(1/2) 



OS 
d\a\ 



\a\[l- 



(E + V + a ) 2 - \a\ 21 







0, 



% OS 



1 



E + V + ib 



. 



(l/2)|£ = -Vo 



1 



a + V + E 



dV ' E + V + ibo (E + V + a ) 2 - |a| 2 

Thanks to ( |VTTg ) & ( jVTTg ), (JVTT§) can be written as 

V + ib -a = 0, 

so that, thanks to ( |V.15| ), 

ib 

Then two cases occur: 



E + a 



Case 1 : \a\ = 0. Substituting in (|V.13|) leads to 

(E + 2a )(a 2 + a E + l) = 



(V.13) 
(V.14) 
(V.15) 

(V.16) 

(V.17) 
(V.18) 
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giving two solutions 



-E±i^A-E 2 , —E ± i V4 - E 2 Tr , 
a = ^ , ^o = \ , V = 0, (V.19) 



oo = = Vo = ^F". (V.20) 

Case 2 : |a| ^ 0. Then, flXTTJ), (JV7LJ) & QVT7D give 

V + E V + E 2 3V + £ 
oo = 2 — ' l a ' = ( — 2 — ^ ~ 2 ' 

3^ 2 - 2£?y + 4 - £ 2 = . 
It leads to the following solution 



E ± 2iy/3 - E 2 
Vo — o 



2E ± zV3 - -E 2 
«o = § 

z6 = —E =f iV3 - £ 2 , 



5E 2 - 12 ± i4EV3 - £ 2 



|a| 2 = ^—^ . (V.21) 



As N — > oo the leading contributions are given by the value of the action at 
the saddle, and we find, for the first saddle point 

S = . 

For the second saddle point, if E 2 < 4, 

S = ^^ + 4ln(\E\/2)<0, 

and, for the third saddle point, if E 2 < 3, 

S = £ 2 /3-ln3±2i(-0 + sin20), < , 
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where E/y/3 = cos0, y/3 — E 2 /y/3 = sin0. Consequently, in the large N 
limit, only the first saddle point contributes to the density of states, and 
gives the usual semi-circle law. Therefore our model as expected is in the 
class of universality of the semi-circle law as iV — ► oo. The saddle points 
are degenerate and one should bound the difference between saddle point 
contribution and the correction corresponding to out-of-the-saddle part of 
the integral. This should be done e.g. like in pO . 



VI The diagonal quartic terms 

In this section, the corrections due to the B term are shown to be small as 



N — > oo. As noticed in Section [TV], the bosonic part of the quartic term B has 
the wrong sign. For this reason before integrating over superfields it will be 
convenient to write a Taylor expansion with integral remainder successively 
for each of the N sectors appearing in the sum for B: 

B = J2 B « ■ 

a 

To first order the expansion gives 

e~ Ba = 1 - I B n e~ Ba e +tBa dt . 



o 



This Taylor expansion either suppresses B a from the exponential of the action 
or generates a remainder term R a — — f Q B a e~ Ba e +tBa dt. Let p(N) be the 
integer part of N/ y/\ogN. The expansion is stopped to the order p(N) = p. 
This gives: 



e- B 



1+ Rp > 



PC[l,-.,lf] 
0<\P\<p(N) 



where 



R P = n r <* if i p i < p( n ) > 

R P =Y[R a \ J e- Ba if \P\ =p(N) . (VI.22) 



a&P o>maxP 
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Let Q be the complement of P in [1, N\. The term 1 was treated in the 
previous sections. The remainders terms Rp must be shown to be 0(1/N) as 
N — ► oo. Let Rp(N) = T2pc[i n] Rp be considered first. For this term it is 

\P\=p(.N) 

not necessary to perform any saddle point analysis and it is sufficient to return 
to the original expression given in (|IV.4|) , hence to a functional integral over 



V. However some e~ Ba -factors are missing or appear with reduced weights. 
To correct for this, all the e~ Ba -terms are recombined with the A-term to 
reproduce the initial functional integral (|1V.4|). However, there are quartic 



correction terms e tBa or e Ba . The important remark is that the bosonic part 
of these terms has now the right sign ! Therefore they can be represented as 
a well defined integral over a new auxiliary field W a . For instance 



e -2tS+S- a S+ a S a _ e -2t\S a S- a \ 2 _ J dW a d\y ae -\ w a\ 2 e tV2iW a S a S- a +cc 

With slightly condensed notations, this leads to 

R p(N) = I S + s(l[ f\ a dt\ m ^ E - v+VTtw ^d^+d^{V,W) . 

FC[1,...,JV] J ^QGP^° ' 

\P\=p(N) 



Then a complex translation Vq h- > Vo±zA _1 is performed, with the same sign 
as the imaginary part of E in order to avoid crossing of singularities. In other 
words 



/+0O P + OO 

e- y °F(V )dV = / e~ v ° e~ 2iX ~ ly o e x ~ 2 F( V + iX^dVo . 
■oo J — oo 

The functional integral can now be bounded by its absolute values every- 
where, namely the following contributions are bounded 

- by 2 N , for the sum over P, that is the total number of subsets of [1, N]; 

- by 1, for the integrals such as J dt; 

- by 1, for the oscillating imaginary integrals; 

- by Gram's inequality for fermions or the Schwarz inequality for bosons, 
for the remainders terms B a ; 

- by 1, for every propagator since the imaginary translation in Vq has 
created an imaginary part proportional to the identity in the denominator of 
the Green's function. 
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This means that each B a term give rise to a small factor A 2 = 1/N for each 
sector in P, hence a factor l/N p ( N >. The two source terms are bounded by 1. 
The normalization determinants are then easily bounded by c N , even without 
using the supersymmetry cancellations, since the operators considered are 
bounded in a finite 2N dimensional space thanks to the imaginary part of E 
which is no longer infinitesimal. Combining all factors leads to 

showing that this correction term is small indeed. 

It remains to prove that the sum of the terms Rp with 1 < \P\ < p(N) 
is small. Let p — \P\. To bound these terms a mean-field analysis will be 
performed (see Section [IV] & |VJ) in terms of the a and b fields, but only 
for the sectors of the theory in the complement Q of P. The functional 
integral to be bounded for a single term is ( |V1.22| ). Now, the e~ Ba, -terms 
are recombined only for a G P with the A term to reproduce the initial 
functional integrals over the V fields and the terms e tt3a , also for a G P, are 
again given by defined integrals over new auxiliary fields W a . Finally, the 
quartic terms, with sector sums reduced to Q, are treated exactly as in the 
previous section, hence a mean-fields a, a , b are correspondingly introduced. 
This leads to a representation 

Rp = / dVp^p dVp t Q d 2 a dVo da dbo Y\ Ra & Cq , 

where Q is the complement of P in [1, ...,N]. In addition, Vpp is the part 
of the matrix V corresponding to rows and columns in P, including the 
new fields of the W type. Vp t Q correspond to one entry in P and the other 
in Q, and the mean field computation is now retricted to Q. The integral 
over superfields gives rise again to a superdeterminant and an additional 
corrections, of the type 

exp [Trlog(l + CV)] , 

where C = R~ Y is the matrix for the N — p analogous problem, as in the 
previous Section |TVJ (see eq. ( |1V.1(J| )), and V = Vpp+Vp^Q is the perturbation. 
This correction term is bounded by 

| exp [Trlog(l + K)]\ < exp (TrK + K* + KK*) . 
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Evaluating C = R' 1 at the saddle point costs a factor exp (Np/2N) = e 2p 
at most. Each term B a naively gives a factor A 2 = 1/(2N) when evaluated, 
but this simply compensates the sum over a when p is nonzero but small, 
so we have to gain an additional 1/N factor. Adding a few expansion steps 
gives such a small additional factor 1/N (already for the first non zero value 
p — 1). This can be seen by integrating by parts the superfields in the vertex 
B a which has been taken down the exponential. By supersymmetry, the 
vaccum graph corresponding to a contraction of the four fields at the vertex 
vanishes. This is absolutely necessary since this graph by simple scaling is 
proportional 1/N and cannot have any additional 1/N factor. Its vanishing 
can be checked by hand: 

- the self-contractions of the bosonic piece 2S+S- a S^ a S a give a factor 



give a factor -2 (since there is one fermionic loop giving the minus sign and 
one sum over a giving a factor 2 only after the contractions); 

- the selfcontractions of the term | J2a(^ta^^a) 2 ar e clearly supersym- 
metric and also add up to 0. 

Consequently, at least one field of the vertex B a has to contract to the expo- 
nential. Performing two contractions in turn, gives always at least a factor 
p/N 2 at the end instead of the naive 1/N factor. Indeed the worst case 
corresponds to the non trivial contraction term being of the type Vpq. This 
generates a new factor 1/N but a new sum over (3 G Q, which costs N—p ~ iV 
so nothing is gained yet. But contracting this (3 field again either generates 
a diagonal term, hence a new 1/N factor, with no new sum, or returns to a 
Vq^p term. This last situation generates a new 1/N factor and a new sum 
over sectors 7 but this time this new sum is restricted to P, so it costs only 
a factor p instead of N\ Hence at worst, after these two contraction steps, a 
total factor p/N 2 instead of the naive factor 1/N is associated to each vertex, 
as announced. Now the sum over P costs a total factor N\/p\(N — p)\, hence 
is bounded by c p [N/p] p . Combining all factors, the sum of contributions of 
such terms with p 7^ is bounded by ^2p=i[c/N] p < c'/N. It is therefore at 
least as small as 1/N. 

This technique could also apply to the term C in the formula ([IV. | ) . Indeed 



+2; 



- the self-contractions of the boson-fermion piece 




a' 
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its presence does not modify the large N density of states. This remark 
might be used to simplify section 3. However it is quite natural to consider 
the diagonal part of V on the same footing that the other mean field terms 
in the supermatrix R. 
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